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1 Introduction 



The study of semileptonic and leptonic decays allows the determination of funda- 
mental parameters of the Standard Model (SM), which are related to the elements of 
the Cabibbo-Kabayashi-Maskawa matrix. Examples of these processes are shown in 

Fig.m 
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Figure 1: Diagrams for the reactions — > /U+z/„ and D' 



K n+u^. 



The strong interaction of the participating quarks introduces a complication be- 
cause its effect is challenging to estimate accurately. Theory's difficulty in calculating 
phenomena rooted in strong physics seriously hamper some measurements related to 
weak physics in the B sector that need such input, for example extracting Vub from 
semileptonic and leptonic B decays, or Vtd from B^ — B^ mixing. Theory tools have 
become available, but they require calibration or verification. The similarity of the 
charm and the bottom quark allow to test the same methods in D decay, where the 
quantities in question can be determined to excellent accuracy because the corre- 
sponding CKM matrix elements are well-constrained from other sources. In addition 
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to being of interest in their own right, experimental charm results thus have a much 
larger impact on the field. 

This presentation discusses Icptonic and semileptonic D and Dg decays, where 
results from the BaBar, Belle, CLEO, and FOCUS collaborations are reported. The 
challenge in all these cases, aside from signal purity and background concerns in 
general, is to cleanly identify decays with neutrinos and to determine the kinematic 
properties of the undetected particle through external constraints. 

The CLEO experiment benefits from the clean experimental environment that 
arises from running at or slightly above production threshold, e^e~ ■?/'(3770) 
DD or e+e~ — > D*Ds — > jDsDg. Identifying one of the D^gj by reconstructing its 
decay to a well- identified final state ("tagging") already guarantees a second D[g] in 
the event, which is then analyzed for a particular reaction. Because of the known 
final state in the first step, it is also possible to not require a tag, but to add all 
reconstructed tracks and showers to infer the neutrino momentum. Finally, and of 
particular importance, the experimental set-up also allows to count the number of 
produced decays accurately, thereby facilitating absolute normalization. 

The B factories usually use the continuum process e'^e" — > cc at center-of-mass 
energies near the T(4jS'). They obtain D mesons via the decay chain D* — * t^sD, 
where tt^ is a slow pion to mark the D flavor. It is possible to analyze both 
and in this way, but in practice, this method is mostly applied to D*~^ '^t-^^i 
which also has the most favorable branching fraction. Since determining the number 
of D decays produced in fragmentation is difficult, measurements are often performed 
relative to another D^^] decay mode. Tagging in this case means that the presence of 
one reconstructed D meson from e^e~ — cc does imply a second charm particle in 
the event (which then is the signal side) . Even in the absence of a way to normalize 
the rate, it is still possible to compare the shape of kinematic distributions with 
theoretical predictions, for example form factor shapes (see below). 

2 Leptonic Decays 

Leptonic D decays proceed through annihilation of the constituent quarks into a W ^ 
followed by its decay into a lepton and the corresponding neutrino. The partial width 
for i = e, 11, and r is given by 

r{D-^ tu,) = ^GlflmjMn (l - ^) Vc/, (1) 

where fu is the D meson decay constant and and Mjj are the lepton and D meson 
masses. Gp is the Fermi coupling constant. For Dg decay, the mass of the Dg, 
the decay constant of fos, and Vcs are used instead. A similar formula applies for 
B^gj decays. The quantities determined in leptonic D^g^ decays are the charm decay 
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constant, given that the masses and the CKM elements are well-measured. The 
decay constant describes the strong interaction surrounding the annihilation process. 
Experimental accuracy in the few-percent range is crucial to validate calculations of 
the decay constant in both the charm and the bottom sector. Further goals are: 
comparing the lepton species to put limits on physics beyond the Standard Model, 
and determining fos/fDi where some technical uncertainties cancel. The Standard 
Model predicts the following ratios for the production rate of r : yU : e: 

D ^ iur. 9.72 : 1 : 0.00002, ii^e : 2.65 : 1 : 0.00002. (2) 

In principle, it is possible to determine the decay constants from decays to all 
three lepton species. A few experimental considerations to take into account: — > 
e'^i'e is beyond current experimental reach (low branching fraction due to helicity 
suppression). Dj^j — » r+i/^ is most copiously produced, but the observed rate is 
lowered, and the measurement complicated by, the subsequent decay of the r, which 
involves at least one more neutrino. D^^^ P-^i^fi is therefore the experimentally 
accessible. 

2.1 Leptonic D decays 

CLEO studied the decay /^^^^ in 281 pb~^ of 1^(3770) data using the tagging 

technique [Ij. The decay is identified with a single muon-like track on the signal 
side, and candidate events are required to have a missing mass squared near zero. 
The missing mass squared is calculated using the beam energy -Ebeam, the muon 
energy and momentum Ef^+ and and the momentum po of the D: MM"^ = 
) - {-PD- -Pn+) , which for signal events corresponds to the neutrino 
mass. The signal distribution is shown in Fig. [2J Fifty signal candidates are found, 
with 2.8 background events (mostly decay to tt'^tt^, r^(— > vr"*"z/), and K^n'^ tails 
from the well-separated peak at higher MM^) expected. This allows measurement 
of the branching fraction B{D^ fJ'^^fi) = (4.4 ± 0.7 ± 0.1) x 10^^ and hence 
fn = (222.6 ± 16.7;3;4) MeV, with Vcd as an external input. With the SM ratios 
from Eqn. [2] applied to the measured branching fraction B{D^ —>■ yU"*"//^), one expects 
B{D+ e+z/e) ~ 1 X 10"^ and B{D+ r+z/^) ~ 2 x 10"^ Slight modification 
of the signal side selection criteria gives access to i = e, r, by asking for either an 
electron-like track (no events seen, B{D^ —>■ e+z/g) < 2.4 x 10~^), or by asking for a 
pion track and a missing mass shifted away from zero, signalling the decay t ^ nu 
(no significant signal, B{D~^ —>■ r^Vr) < 3.1 x 10~^). 

2.2 Leptonic Dg decays 

Both BaBar and CLEO have recently studied leptonic Dg decays. BaBar bases their 
study on 230 fb"^ of data in the T(45) region. They use kinematic constraints, require 
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Figure 2: ^'^i'^'- reconstructed neutrino mass for signal events [T]. 



an identified D+, or D+, and use the signal side decay chain D*+ — >• 

H'^v^ to provide a clean Dg sample ^ (extra particles may be present in 
the event). The neutrino momentum is inferred from all other measured reaction 
products. The signal distribution is the D* — Ds mass difference m(7/i+z/^) — m(/i+z/^), 
shown in Fig.[3l Since the D* production rate is not precisely known, BaBar determine 
the ratio relative to {K^ K~)tt^ with m{K^K~) in the mass region {±2r^), 

r(D+ ^ /i+i/^)/r(D+ ^ 07r+). With B{D+ 07r+) = (4.71 ±0.46)% ^ they arrive 
at B{Df /i+i/^) = (6.74 ± 0.83 ± 0.26 ± 0.66) x lO^^, where the third error is due 
to the normalization uncertainty. This leads to fr,^ = (283 ± 17 ± 7 ± 14) MeV. 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 
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Figure 3: BaBar's signal distribution AM = m{'-)^^v^) — m{^'^v^) for the decay 
D+ /i+z/^ detected in D*^ -fD^ [4J. The structure at AM ~ 0.07MeV/c2 is due 
to the decay D* t^^Ds] the peak underneath the signal in the background shape is 
due to TUr- 
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CLEO's analysis relies on ~ 200 pb of data taken at a center-of-mass energy of 
4170 MeV. The decay chain is e+e" D*Ds {jDs)Ds. One Ds is tagged, the 
other is analyzed for the signal decay. Two strategies are used. 

A. As above, the missing mass squared is examined to identify the signal. The 
lone signal side track is required to not be an electron but instead to be muonic or 
pionic. An energy requirement of 300 MeV accepts 99% of all muon tracks and 60% of 
all pion tracks. We distinguish cases with an energy deposition below 300 MeV in the 
calorimeter (typical for f^~^i^n, but possible for r —>■ vrz/^) from those that have 
above 300 MeV (in which case the track is likely to be a pion from r — ttUt-). The data 
show a clear enhancement in the expected D'^ /^^^'a* signal region. The Z)^ — > rz/^ 
events are more spread out due to the presence of the additional, unreconstructed, 
neutrino, but signal events are seen as well. The signatures for i = /i and r overlap, 
hence the summed distribution is fit for the two branching fractions, where the SM 
ratio for the two is assumed. CLEO's preliminary result, combining the two, is 
B{D+ /i+i/^) = (0.664 ± 0.076 ± 0.028), or fo, = (282 ± 16 ± 7) MeV. If the track 
is required to be consistent with an electron, no candidates are found, resulting in an 
upper limit of B{Df e+u^) < 3.1 x 10"^ 

B. The second approach, based on the same data, uses the decay chain r —>■ eV^.v^. 
Its product with the — > t^Vj. branching fraction (~ 6 — 7%) is about 1.3%, to be 
compared with the inclusive semileptonic branching ratio D"^ —>■ Xe+z/g ~ 8%. The 
analysis procedure demands a sole electron-like track on the signal side, and limits the 
energy not associated with the other identified decay products in the calorimeter to 
be less than 400 MeV. No additional energy deposition is expected for signal events 
other than the transition photon from D*Ds {'yDs)Ds, upon which no selection 
requirements are placed, and showers resulting from interactions of the decay products 
of the tag side with the detector material. The signal distribution, together with 
background estimates from Monte Carlo simulations, is presented in Fig. HI This 
analysis leads to B{Df ^ r+Ur) = (6.3 ± 0.8 ± 0.5)% and /d, = (278 ±17 ±12) MeV 
(both preliminary). Since the two measurements are complementary, one can form 
the average = (280±12±6) MeV and also use them to measure the ratio B{D^ — >• 
r+i/^) : B{D^ /i+z/^) = (9.9 ± 1.9), consistent with the SM expectation of 9.72. 
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Figure 4: CLEO data on leptonic Dg decays (preliminary). Left: Missing-mass- 
squared distributions from data corresponding to (i) D'^ fJ'^i^^i + t'^i^t, (ii) — ^ 
t+i/t-, (iii) — >■ eUf,. Right: Df t^Vt-: Energy deposited in the calorimeter for 
tagged Dg events and a single electron-like track on the signal side not accounted 
for by the tag or the electron candidate. The circles are data; the curves are MC 
predictions for signal as well as several semileptonic background sources. 



2.3 Leptonic Decays: Summary 

A visual comparison of current experimental and theoretical progress on the decay 
constants is given in Fig. [51 In summary: /^j from — >■ /i^t'^ is measured to a 
total relative error of 8%; the decays to e and r are presently not experimentally 
accessible, although the current upper limit gives rise to the hope that a signal will 
soon be within reach for r. The decay constant fr,^ has been measured in — > /U+z/^ 
(BaBar, CLEO) and r"*"//^ (CLEO) to an accuracy of 5-8%. The precision of 

theoretical calculations are in this region as well. However, the dominant experimental 
uncertainty is still statistical, which is - in principle - much easier to improve upon 
than systematic errors such as those theory is in the process of overcoming. 
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Figure 5: Recent experimental and theoretical results on Jd and Jd^. 



3 Semileptonic Decays 

Semileptonic decays D —>■ hiu^ provide yet another scenario within which to study the 
impact of strong force. The underlying weak process at the quark level, for instance 
c — > W*q with W* — > eu, can be calculated, but the observed rate is modified by 
the QCD interaction between the participant as well as the spectator quarks, which 
clouds the simple picture. 

3.1 Branching fractions 

The current experimental uncertainty on the D ^ K, vrez/ branching fractions is 
in the percent regime, thereby posing a challenge to the precision of the LQCD 
predictions [12]. 

The search is on for rarer modes, and CLEO have improved considerably upon 
previously achieved accuracy in many modes (Tabled]). — > rje'^u has been ob- 



7 



served for the first time (> 5cr), as has the first multi-body semileptonic decay: 
K^-K^Ti^ e^Ve (> 4cr), which is found to mostly proceed via ^1(1270)^ K^7i^7i~. 
The uncertainties on uje^u branching fraction have been halved, and the upper 

limits on the branching fractions — »• rj', (pe^u have been tightened by two orders 
of magnitude. Once an rj'ei'e signal is observed this will allow a comparison with 
predictions for the ratio D ^ r] : D —>■ rj' . 
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T: 8.86 ± 0.17 ± 0.20 8.9 ± 0.4 6.8 ± 0.8 
U: 8.75 ±0.13 ±0.30 

T: 3.58 ±0.05 ±0.05 3.41 ± 0.09 3.43 ± 0.14 
U: 3.56 ±0.03 ±0.11 
Belle: 3.45 ± 0.07 ± 0.20 

T: 0.397 ± 0.027 ± 0.028 0.44 ± 0.07 0.31 ± 0.15 
U: 0.383 ±0.025 ±0.016 

T: 0.309 ± 0.012 ± 0.006 0.27 ± 0.02 0.36 ± 0.06 
U: 0.301 ±0.011 ±0.010 
Belle: 0.255 ± 0.019 ± 0.016 

T: 0.232 ± 0.020 ± 0.012 0.24 ± 0.04 0.31 ± 0.06 

T: 0.156 ±0.016 ±0.009 0.19 ± 0.04 
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Table 1: CLEO's preliminary [6] and Belle's published [15J measurements for 
D semileptonic branching fractions, and comparison with PDG06 ^ as well as 
PDG04 [8]. The PDG06 results include CLEO's pubhshed numbers based on 56pb~^ 
of data taken at the ?/'(3770). "T" and "U" stand for CLEO's "tagged" and "un- 
tagged" analyses, respectively, and the corresponding entries are not to be averaged 
because of sample overlap. CLEO's results are all obtained with £ = e. Belle's with 
i = e or p, and the PDG numbers are averaged over e and p, where available. 
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Figure 6: Preliminary CLEO results: Signal distributions for — > rje'^i' for rj ^ 'j'j 
(left) and r] tt+tt^tt^ (middle); D° — K^Tr+vr^e+z/. 



3.2 Form factors 

3.2.1 Decay to a pseudoscalar 

The dynamics of semileptonic D —>■ hiu^ decays can for the case of a pseudoscalar 
hadron in the final state in the limit of small lepton masses be described as follows: 



fliq') xplx\V,,\', (3) 



where is the momentum transfer to the W*, f^{q^) is the form factor function 
describing the probability to end up with a hadron of type h in the final state for 
a given g^, ph is the momentum of the outgoing hadron, and Vcq is the appropriate 
CKM matrix element. 

Theoretical predictions for the form factor shape can be tested against the dis- 
tribution in data, corrected for the g^-dependent detection efficiency. This is detailed 
further below. The normalization is determined by the product \Vcq\ x /+(0). Due to 
the precision with which Vcs and Vcd are determined in other experiments (2 — 4%), 
a measurement of \Vcq\ x /+(0) determines /+(0). 

Examples of form factor determinations are displayed in Fig. [71 The simplest 
reasonable parametrization of the pseudoscalar form factor is the single pole shape, 
~ 1/(1 — g^/MpQig) [To], where the nominal setting of the parameter is Mpoie = 
M[Dl)[M{D*)] for h = K[7i]. More sophisticated models are the modified pole 
model, 

^^^^^ (l-«gVMp^ole)x(l-gVMp^ole)' ^ ^ 
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or, recently, the Hill series parametrization [TTj. BaBar [13], Belle [15], CLEO [6j, 
and FOCUS [H] compare their various results for D ^ n, Kiu data for D = D^'^ 
with the models and find reasonable agreement with all of them. In particular the 
Belle and CLEO measurements boast superb statistics and an excellent resolution. 
In general it is found that unquenched LQCD predicts the form factor shape a little 
too high especially for D — > vr [6l |15] , but has achieved useful uncertainty, as visually 
evident from the agreement between the data points and the curve obtained from an 
interpolation between LQCD predictions at several points. Fig. [3 Fitting the pre- 
dicted LQCD data points to the modified pole model shape, thereby extracting /+(0) 
and a, allows a comparison within this model between theory and the experiments. 
This is presented for /+(0) in Fig. [HI 




q^(GeV^) 

Figure 7: Form factor shapes. Left (from Ref. BaBar {D —* KeUe, points) and 

FOCUS ^14] (triangles, D Knu^) data and an unquenched LQCD calculation 
(band). Right (from Ref. [U]): Belle results on D ^ K (top) and D ^ n (bottom), 
with unquenched [12] and quenched [TB] LQCD predictions overlaid. 

Alternatively, using f^^^ = 0.64 ± 0.03 ± 0.06 and f^^^ = 0.73 ± 0.03 ± 0.07 
from LQCD [12], the CKM matrix elements are found to be in good agreement with 
current world averages. However, the uncertainties (~ 10% relative) on the form 
factor magnitudes dominate. This constitutes a check of the LQCD calculation. 

3.2.2 Decay to a vector 

For a vector hadron in the final state, such as D ^ K*ev (Cabibbo-favored) or Z) — >■ 
pep (Cabibbo-suppressed), more form factor functions enter the stage. The decay 
amplitude can be described in a parameter-free way using helicity basis form factors 
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Figure 8: Comparison of form factor normalization for D ^ K and D ^ tt obtained 
by unquenched LQCD [1^, two CLEO [6| analyses, and Belle [T5] . 



as specified in Ref. [T7j. However, the helicity form factors are hard to calculate; a 
more traditional approach is to just assume spectroscopic pole dominance and cast 
the expression for the amplitude into linear combinations of the following functions: 

M^') = , "^fL i^ = l,2),T/(g^) = r^^- (5) 

1 - q^Mi^ 1 - gVM^ 

The pole masses are often fixed to be My = 2.1 GeV, M^^ ^ = 2.5 GeV. If one then 
defines Ry = V{0)/Ai{0) and R2 = A2{0)/Ai{0), one ends up with only two free 
parameters, albeit at the expense of an assumed shape. The parameter-free approach 
in Ref. [17] circumvents this. 

An analysis of K'lx^e^v CLEO data uses a projective weighting technique, 

by which the expected shapes from helicity form factor contributions are fit to the 
data, thereby extracting the (bin-wise) amplitude for each of the form factors [18j. 
When the spectroscopic form factors as given in Eqn. [5l translated into helicity form 
factors, are overlaid, good agreement is found (Fig. [9]). Additional conclusions are 
that a term describing interference with a non-resonant s-wave Kt^ component is 
necessary, and that no evidence is found for d- or /-wave contributions. 

A recent preliminary study by BaBar [19J of D"^ (pe'^u, using about 13 x 10^ 
signal events in 78.5 pb~^, not only exhibits a beautifully precise measurement but also 
confirms an earlier FOCUS result. This actually resolves a controversy: The ratios 
Rv and R2 for K*e^u and (pe^v are expected to be similar because the 

CKM matrix element involved is Vcs in both cases, and the only remaining difference 
is then the spectator quark flavor {d vs. s). Agreement within 10% is expected. This 
was not borne out by earlier data, which was not inconsistent for Ry but disagreed 
for R2, except for a FOCUS measurement |20] that had R2 in agreement with the 
world average of -Rf^*. 
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Figure 9: D ^ Kireu [18]: Four parameter-free helicity-basis form factor products as 
extracted in CLEO data, overlaid with form factors constructed from spectroscopic 
pole model shapes using Rq = 1.505, R2 = 0.875, and s-wave parameters A = 0.33 
and 5 = 39°. 



CLEO has presented preliminary results for the first Cabibbo-suppressed pseu- 
doscalar vector form factor measurement, D peu. A clear signal for both isospin 
states is observed, making it possible to bin in kinematic distributions. The branching 
fractions from the 281 pb~^ data sample are found to be in good agreement with pre- 
vious results (which are dominated by the CLEO-c results from the 56 pb~^ sample, 
see Tabled, as is the partial width p'e+u = (0.41 ± 0.03 ± 0.02) x 10"^. A 

simulatenous fit to — >■ p^e^Ve and p^e^Ve (linked through isospin) results 

in the following determinations: R^r = 1.40 ± 0.25 ± 0.03, R^ = 0.57 ± 0.18 ± 0.06. 

D —>■ pev and D — > K*ev studies are important to the B sector for example by 
aiding in the extraction of iKblVl'^i^P through dV{B pep) : dV{B K*U) [9J, 
which requires either measurement of the B ^ p and B K* form factors or a 
calculation validated by the corresponding D decays. 

3.3 Inclusive semileptonic decays 

CLEO has determined the semileptonic inclusive branching fraction D — > Xe'^Ue for 
both charged and neutral D mesons [21]. The measured branching fractions B{D^ — > 
Xe+Ue) = (6.46±0.17±0.13)%, B{D+ Xe+^e) = (16.13±0.20±0.33)% agree well 
with the sum of all exclusive modes, although there is room for as-yet-unobserved 
exclusive decays at the level of i3 ~ 10~^. Isospin symmetry is observed within errors 
for the inclusive semileptonic partial widths, F^i+ZF/jo = 0.985 ± 0.028 ± 0.015. 

Aside from the branching fraction, a quantity of interest is the electron momentum 
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Figure 10: Comparison of form factor parameters R2 (left) and Ry (right) for D 
K* and Dg ^ (p semileptonic decays, from Ref. [19]. The data points all represent 
Dg ^ (j) analyses. The lines give the ilex band for an average of experimental 
D ^ K* results. 



spectrum. Without weak annihilation contributions {D^ — > gW^ Xe^v) agree- 
ment between the spectra from and is expected. Weak annihilation would 
modify the upper end of the electron spectrum. Within uncertainties, the and 
distributions agree (Fig. [TT]) . 
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Figure 11: Electron momentum spectrum in D ^ Xe^Ve for charged and neutral 
D mesons |21j . 
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3.4 Semileptonic decays: Summary 



Many new results on D branching fractions and D ^ P,V and Dg ^ V form factors 
have become available. Of particular interest are the shape and the normalization of 
the form factor functions. Experimental accuracy is at present still consistent with 
most calculations on the market; this may change as the experimental and theoretical 
uncertainties decrease. 



4 Conclusions 

Precise results in the D sector have improved our comprehension of the QCD effects 
accompanying weak interactions and allowed to sharpen theoretical tools. Thanks 
to the similarity of the heavy quarks c and b, common calculation techniques can 
be applied to the estimation of D and B decay properties. Further progress in the 
D sector and consequently the B sector is in sight as data samples with modern 
detectors are being enlarged. 
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